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ABSTRACT 


Oceanographic and meteorological data from Ocean Weather 
Station (OWS) NOVEMBER were analyzed by statistical and 
spectral methods in order to describe the nature and perio- 
dicities of the air-ocean energy exchange. Salinity data 
from 1968 through 1970 were compared with daily observed 
parameters commonly associated with changes in the salinity. 
Regression analysis showed surface salinity to be highly 
correlated with surface pressure which is interpreted as the 
surface salinity responding to a baroclinic transport. At 
this location none of the classical parameters usually 
associated with fluctuations in the surface salinity were 
highly correlated. 

Twenty-four years (1947-1970) of daily averaged, 
observed climatological values, as well as heat exchange 
terms computed from them, were compared with the sea-surface 
temperature time saisiec and spectrally analyzed. The 
Spectral analysis yielded spectral density, phase and 
coherence functions indicating the response of surface 
temperature to the heat exchange terms. From these it was 
evident that at the yearly cycle the net heat exchange, which 
is shown to be equal to the horizontal heat advection, is 
of the same order of importance as solar radiation in its 


effect on the fluctuation of sea-surface temperature. 
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I. INTRODUCTION 


A. OBE Ue 

Sea-surface temperature is one of the most important 
parameters in any study of air-sea interaction. Anomalies 
of this parameter, measured over periods of a week to 
several months, have been linked to gross atmospheric 
effects. However, few studies have investigated the time 
response of the sea-surface temperature over long periods 
to those factors which govern its fluctuations. 

Until the more recent decades such a study would have 
been limited by both instrumentation and general paucity 
of data. This resulted in reliance upon monthly or annual 
means to describe the state of the ocean. The advent of the 
ocean weather station program in the 1940's, and recently, 
the use of satellites has greatly increased the amount and 
frequency of data available to the researcher. Also, the 
advent of high speed computers has made the use of time 
series analysis routines relatively easy and inexpensive to 
use. 

The past limitations made statistical studies of changes 
of surface temperature on time scales of a day or less, 
difficult. However, daily fluctuations can be quite inform- 
ative. If advection is ignored, the sea-surface temperature 
at OWS NOVEMBER was found to fluctuate primarily in response 


to the day to day changes in air temperature, the amount of 





insolation, wind-wave mixing, and the amount of cloud cover 
as well as the elevation of the sun. The magnitude of 
these diurnal fluctuations can be as large as the monthly 
fluctuations. 

Sea-surface temperature fluctuations are important in 
many oceanographic parameters, such as the stability of the 
water column. This, in turn, is of great consequence in 
sound propagation. Anomalies of sea-surface temperature, 
likewise, have been linked to climatic variations over the 
ocean. 

The salinity is important in that it affects the 
density of sea water. Thus, it affects both vertical 
Stability and sound propagation. Together with the water 
temperature, one can infer transport characteristics of a 
region based on water mass analysis. 

This study examines changes occurring at the surface 
of the ocean over time scales varying from one day to 
several years. The two parameters of primary interest are 
the sea-surface temperature and surface salinity. The area 
of study is Ocean Weather Station (OWS) NOVEMBER located at 
30° North and 140° West. The data used are synoptic data 
morected by the United States Coast Guard at this position. 
The data are of two types: hydrographic Nansen cast data and 
shipboard meteorological observations. 

The thesis is organized to first present correlations 
of the primary parameters versus various meteorological 


parameters. These correlations provide an indication of 


(as: 





those parameters which can be used as predictors for the 
sea-surface temperature and salinity. 

Secondly, a spectral analysis of the data is made. 
From this, the synoptic atmospheric mechanisms which cause 
the fluctuations of the sea-surface temperature are identi- 
fied. Periods of high spectral energy, along with their 
corresponding phase and coherence relationships are examined 
to gain insight into the time scales associated with each 


mechanism. 


S. Move WOR eRIE EITERATURE 
Iz The Subtropic Region 

Ocean Weather Station NOVEMBER, located at 30° 
North and 140° West, lies within what Tully (1964) defines 
as the Pacific Subtropical Region. Figure 1] shows its 
location in comparison to other major features of the North 
Pacific Ocean. The region is characterized by three dominant 
conditions: an excess of evaporation over precipitation 
throughout the year, the occurrence of convective mixing, 
and seasonal heating and cooling patterns which result ina 
pronounced seasonal thermocline. 

The excess of evaporation over precipitation (E-P), 
has a distinct effect on the other two traits. Tully (1964) 
estimated E-P to be approximately 80 cm”! yr in the vicinity 
of OWS NOVEMBER. This leads to a removal of fresh water 
from the surface and acts as a driving force for the con- 
vective mixing. This mixing then contributes to the erosion 


of the seasonal thermocline. 
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rae surface Sieur ae on 

The subtropical region of the Pacific is dominated 
by the Pacific Gyre. Figures 2 and 3 show graphically how 
the gyre varies seasonally both in location and intensity. 
In the winter the gyre's axis is displaced northward; in 
summer it is shifted approximately 100 n mi southward and 
intensifies by approximately 0.2 kt. Throughout the year 
the current speed ranges from 0.1 to 0.3 kt. Hansen (1973) 
hypothesized that this intensification occasionally results 
in the advection of large volumes of cold subarctic water 
into the area surrounding OWS NOVEMBER, specifically during 
the period of 1968 through 1970. 

oS The Subarctic Transition Zone 

The North Pacific Gyre, linked witn tne Subarctic 
Gyre, 1S also involved in the formation of a transition zone 
which bounds the subtropical region to the north. Roden 
(1970) first described the subarctic transition zone as an 
area characterized by large temperature and salinity gradients 
and low hydrostatic stabilities. This oceanic front defines 
the boundary between the subtropical and subarctic water 
masses, varying in location between 30° and 40° North. 
Subsequent articles by Roden (1971, 1972) further described 
this zone in the North Pacific. 

The transition zone is largely maintained by the 
Surface wind stress. Its northern boundary is characterized 


by the subarctic water mass, resulting in temperature 
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Summer surface circulation patterns for the North 
Pacific Ocean (Department of Commerce, 1961). 
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inversions, gravitational instabilities, a complicated sound 
velocity structure, and a disappearance of the subarctic 
halocline. It appears to lie in the area of 40N to 45N 
latitude throughout the year. 

The southern boundary lies in the area of 30N to 
32N, depending upon the season, and is displaced to the 
north more in the winter months. It is characterized by 
Sharp temperature gradients in the thermocline and vertical 
instabilities in the mixed, upper layers. A salinity minimum 
1S a common occurrence at this boundary. Also, around 30N 
there is a transition toward the higher salinities and 
higher temperatures of the subtropical water mass with 
Salinities greater than 35°/.. and a strong halocline located 
at a depth between 50 m and 100 m. 

Thus, the transition zone is defined by a deep 
Salinity minimum, at approximately 500 m and a concurrent 
sound velocity minimum of less than 1480 m sec’!. There is 
apparently both a seasonal and yearly variation in the posi- 
tion of the boundary (Seckel, 1968). Both the salinity 
minimum and sound velocity minimum rise as one goes further 
to the north. Roden (1970) also examined both barotropic 
and baroclinic motion in the transition zone and found that 
they rarely exceeded 5 cm cue 
As the transition zone approaches North America, 


it turns to the southeast and becomes associated with the 


heat and salt flux divergence zones noted by Roden (1971) 


We 





in the area of OWS NOVEMBER. These are related to flux 

divergencies which are linked to the presence of horizontal 

advection of salinity as has been noted by Hansen (1973). 
4. Average Annual Heat Balance 

Based upon an examination of the annual heat 
balance of the North Pacific Ocean and the prevalent circu- 
lation patterns, Wyrtki (1965) determined that there is a 
small amount of heat gain in the vicinity of OWS NOVEMBER. 
On the average, he estimated this gain to be approximately 
25 ly day”! (where a langley, ly, = 1 cal cm7*). 

Seckel (1970) computed the heat balance for the 
area Of OWS NOVEMBER during the Trade Wind Zone Pilot Study. 
His studies showed that the average net radiation computed 
on a daily basis was approximately 3 ly day” !. 

In a more recent study, Dorman (1974) studied 20 
years of the heat balance at OWS NOVEMBER and separated it 
into two decades, 195] through 1960 and 1961 through 1970. 
He found the two decades had distinctly separable heat 
balances, with the first decade having a deficit of 


! and the second a 40 ly day~! deficit. These 


90 ly day 
results were based on three-hourly meteorological observa- 
tions. Dorman ascribed this difference to the change in 
climatic regimes described by Namias (1970, 1971), which 
resulted in displacement of cooler waters in the North 


Central Pacific by warmer water, beginning approximately in 


m6 | . 
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5a Surface and Subsurface Temperature Anomalies 


The sea-surface temperature and its anomalies have 
been widely studied in their relationships to both oceano- 
graphic and meteorological forecasting. Emphasis has also 
been placed on subsurface temperature anomalies. Beland 
(1971) illustrated their considerable importance to the flow 
of energy between the sea and the atmosphere. Bernstein 
(1974) used subsurface temperature measurements in the trade 
wind region northeast of Hawaii to identify mesoscale eddies 
drifting westward. 

Namias (1971) in his study of an anomalous warming 
of the northeastern Pacific from May to June 1968, found the 
anomaly extended to a depth of 122 m. The warming, as a 
whole, was attributed to increased insolation and horizontal 
convergence of surface waters with its associated down- 
welling. The warming was identified with the development 
and maintenance of a Pacific anticyclone in June of 1968. 
The depth of the anomaly was attributed to bulk heating and 
downward advection of the warmer surface waters. 

: Clark (1972) concluded that both advection and 
non-advective processes are vital in determining the nature 
of surface temperature anomalies. He hypothesized that 
advection is the more important process in the winter and 
Spring seasons, while the non-advective processes dominate 


in the summer and fall. 


(aa 





Beland (1971) in his study of the North Pacific, 
observed that about 50% of the near-surface temperature 
anomalies extended to at least 100 m depth. Most of these 
were negative anomalies in the months of March and April, 
apparently due to the distribution of heat throughout the 
upper layer by convective mixing. Positive anomalies were 
more common in December through February and were confined 
to the upper 30 m. Whereas the negative anomalies were of 
a magnitude of 1.0°C or more, the positive anomalies were 
generally less than 1.0°C. 

Finally, Favorite and McClain (1973) reported that 
they could observe an orderly trans-Pacific movement of 
areas of anomalously warm and cold surface waters. These 
features were found to be quasi-permanent and apparently 
move around the North Pacific gyre, approximately 180 degrees 
Out of phase and with periods of 5 to 6 years. These 
anomalies have been shown to have a great influence on the 
weather conditions off the west coasts of Canada and the 
United States, as shown by Namias and Born (1970). 

From the above it can be noted that temperature 
anomalies, both surface and subsurface, and whether temporary 
or periodic, inhabit the region around OWS NOVEMBER, and in 
all probability, affect the location of the subarctic 
transition zone and the heat balance through the process of 


advection. 
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6. Recent Studies of Pacific Air-Sea Interactions 

In recent years there have been three studies 
of the air-sea relationships at or around OWS NOVEMBER by 
different researchers using three different methods. Taken 
as a whole, they provide a logical progression, resulting 
in the methods utilized in this study. The first of these 
studies was carried out by Seckel and Yong in 1970, using 
harmonic analysis to study temperature and salinity fluctua- 
tions at Koko Head, Hawaii, and Christmas Island. Hansen 
(1973) used regression methods to study surface and 
subsurface temperatures and salinities at OWS NOVEMBER. 
Dorman (1974) used a spectral approach on time series of 
meteorological values at OWS NOVEMBER, dividing the time 
series into two ten-year decades. Since each of these 
studies utilized a different approach, each will be described 
in some detail here. 

Seckel and Yong (1970) utilized temperatures and 
Salinities measured at two points in the North Central 
Pacific; Koko Head, Hawaii, beginning in 1956 and running 
through 1969, and Christmas Island, beginning in 1954 and 
continuing through 1969. Their hypothesis was that 
fluctuations in the temperature and salinity time series 
reflected air-sea interaction processes (heat exchange, E-P) 
and oceanographic processes (advection, diffusion, etc.). 
Their approach used a Fourier series to provide a least 
Squares fit to the observed values, which was done for each 


year of the available data. This approach expressed the 


ZG 





temperature and salinity as a function of time. Phase 
angles and coefficients -for each year of salinity and temper- 
ature data from both locations were computed and tabulated. 

The results at Koko Head showed seasonal varia- 
bility in the salinity, but not the temperature series. At 
both locations the annual cycle (n=1) was the dominant cycle. 
Also of note were the 60- and 52-day cycles, ascribed by the 
authors to the presence of large geostrophic eddies. fhe 
Christmas Island results were quite similar except that they 
lacked evidence of the eddy-induced seasonal variations. 

The next major method was a regression analysis 
used by Hansen (1973) in his study of temperature and 
Salinity fluctuations at the surface, 50 m, and 200 m at 
OWS NOVEMBER. His data were for 1968 through 1970 anda 
consisted of daily hydrographic records. Deterministic 
approaches showed that the temperature system was seasonal 
So the data were divided into four seasonal records, 
February to April, etc., in three month segments correspond- 
ing to Spring, Summer, etc. Regressions were then computed 
with the surface salinity, and 50 m and 200 m temperatures 
being regressed against the surface temperatures. 

Most of the results showed little correlation, 
except occasionally during the spring and fall transition 
periods. Only a very poor correlation was found between the 


Surface salinity and the surface temperature. From this, it 
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was concluded that these seasonal correlations were due to 
the advection of water by the shifting of the subarctic 
transition zone. 

The final approach was that utilized by Dorman 
(1974) who computed spectra for surface wind speed, air 
pressure, sea temperature, air temperature and dew-point 
temperature for the 10-year period 1961 through 1970. 
Spectra were also computed for each of the four seasons. 
From this it was found that the energy peaked significantly 
at the annual, diurnal and semi-diurnal periods. Variances 
were largest in the winter season and smallest in the 
Summer season, except for the sea-surface temperature which 
had its largest variance in the spring and the smallest in 
the winter. The maximum spectra! density, for all variables, 
was at the annual period. From the annual period the 
spectral densities followed a decrease to the shorter periods, 
where significant spikes at the diurnal and semi-diurnal 


periods were noted in many of the spectra. 





II. TREATMENT OF THE DATA 


A. SOURCES 

The climatological data used in this study were taken 
from the NWRC, Asheville data file of ocean weather station 
observations. These observations were begun at OWS NOVEMBER 
in 1947 and have continued to the present. The observations 
were taken every three hours, except for some short periods 
when all or only some of the observations were not taken for 
some parameters. The data extracted included the sea-surface 
temperature, air temperature, wind speed, surface pressure, 
cloud amount and the present weather codes for 1947 through 
o7 0. 

The salinity data used in this study came from the NODC 
hydrocast data. Since July 1966, the U.S. Coast Guard has 
conducted a continuing program of collecting Nansen cast 
data at OWS NOVEMBER. The surface salinity data for 1968 


through 1970 was utilized in one segment of this study. 


B. TREATMENT OF THE DATA TIME SERIES 

Money Slirface Salinity values from 1968 though 1970 
were extracted from the NODC data set and any gaps in the 
series were linearly interpolated. The daily salinity values 
Pememtnen included with surface pressure, evaporation, solar 
radiation, heat of evaporation, and wind speed, all taken 
directly from or computed from the Asheville data set. To 


this was added the calculated values of precipitation and 
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the evaporation minus precipitation terms to form a three- 
year time series which was then used in the regression 
analysis segment of this study. 

The Asheville meteorological data were used primarily 
to create the 24-year (1947-1970) time series used in the 
spectral analysis segment of this study. The values of SST, 
air temperature, cloud amount, air-sea temperature difference, 
the saturated vapor pressure of the air, the vapor pressure 
over the water, wind speed, surface pressure and wind speed 
Squared were extracted or computed for every three hours 
and then converted into daily average values and these were 
the values utilized in all the spectral analyses. For the 
cloud amount parameter, only daylight observations were used. 
From these data sets the heat budget terms were computed on 
a daily basis. For the occasional gaps in the data, linear 


interpolations were made. 


eS THE HEAT BUDGET EQUATIONS 

Empirically derived formulas were used to create the 
heat flux data for this study. The heat budget equation 
utilized excludes the effects of heat gained or lost by 


ocean currents: 


HEAT CHANGE = HEAT INCOMING - HEAT LOSSES 
QN = QS - QE - QB - QC 


ey 





where QN is the measured heat change, i.e., the heat gained 
or lost by the water due to the exchange processes. On the 
right-hand side of the equation, QS is the incident solar 
radiation, QB is the effective back radiation, QE is the 
heat of evaporation and QC is the conduction of sensible 
heat. Each of these terms will be described in detail below 
as to their formulation. 

The equation used for QS is an improved version of that 


described by Seckel (1970). The basic form is: 


Q = A, “t A, cOos¢ + By sing + A, cos2¢ + Ba, sin2do (3) 


From the Smithsonian Tables (List, 1951), utilizing a 
transmission coefficient of 0.7, the above terms can be 


computed as follows for latitude 30N: 


gps 22.05 + 674.76 cose = 551.7 (4) 

Eee 297.26 cos (2490) = -178.7 (5) 

By = -+- 6.75 + 224.38 sin 2 = 105.4 (6) 

AD = - 1.32 + 16.10 sin 2(2-45) = - 9.4 ion 

Epes «CO -«YC« 04 + «29.76 cos 2(2-5) = 18.1 (8) 
where: 


2 is the latitude 
t is the time of the year in days 


C is the cloudiness in tenths, ranging from 0 to 1, and 
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¢ = AAG (t-21) (9) 


The cloud factor correction utilized was one for 


cumulus type clouds (Seckel, personal communication): 


ena oaieo = 0.363c- 


An albedo factor of 0.95 was utilized. 

The new equation results in values approximately one- 
fourth to one-third greater than previously estimated. 

The effective back radiation, QB, is a Berlyand 


expression utilized by Seckel (1970): 


QB = 1.14 (10)7’ (273.16 + i (11) 


(0.39 - 0.05 ye.) (1 - 0.6e*) 


Pee ceenOy (273516 + — 


where a is the sea-surface temperature in degrees C. It 
utilizes a cloud coefficient of 0.6 and neglects latitudinal 
variation. 

QE, the heat of evaporation depends upon the difference 
between the saturation vapor pressure over the sea at the 
sea-Surface temperature and the vapor pressure of the air 


above the sea surface: 
QE = 3767 Ch (0.98 erence) W (olic) 
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< is the vapor pressure of the air (mb) 


e is the saturation vapor pressure over water 
at sea-surface temperature (mb) 


D is the drag coefficient 


W is the wind speed (m sec”! ) 


Whereas Seckel utilized a varying drag coefficient, this 
meeay Utilized a constant drag coefficient of 1.3 x ily, 
As in Seckel's study, the vapor pressure difference was 
multiplied by 0.98 to account for the salt effect. 


The final heat balance term is the conduction of 


sensible heat, QC: 
OG 64 7 cole) ah (13) 


where: 
i, is the sea-surface temperature (C) 
T, is the air temperature (C) 


W is the wind speed (m see!) 


It can be seen readily that QC is dependent upon two factors, 
the air-sea temperature difference and the wind speed. It 
is usually the smallest of all the heat exchange terms. As 
with Seckel's study, a conduction factor of 4.0 is used. 

The above equations were those used to create the five 
heat flux time series used in this study. QS was computed 


using only the daylight cloudiness observations, night 
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values being ignored. A daily value for each day was 
calculated for the span of 1947 through 1970. They were 
compared with Seckel's (1970) results to guarantee their 
accuracy, where possible. The utilization of the newer QS 
formula resulted in higher values for this term than he 


reported and correspondingly more positive values of QN. 


D. GOMPUTATION OF PRECIPITATION 

In relation to ene study of the surface salinity at 
OWS NOVEMBER, it was vital to have some measurement of E-P, 
the evaporation-precipitation term. The evaporation was 
readily available, having been calculated from the data 
already available. However, no direct measurement of pre- 
cipitation is made routinely at ocean weather stations. 
To gain a measure of the precipitation, Tucker's (1961) 
method was used. Elliot and Reed (1973) applied this method 
along the Pacific Northwest Coast with good results. This 
approach is based on the current weather code which gives a 
numerical value to the type of weather being observed at 
three hour intervals; this yields an estimation of the amount 
of rainfall. The amount of rainfall was based on the amount 
of precipitation observed for a similar weather code at a 
land station where actual rainfall measurements were avail- 
able. Based on a frequency analysis for all the codes and 
the measured rainfall, Tucker was able to fit a numerical 
System of equations to the prediction of precipitation, 


using only the weather code as an input. 
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In Tucker's system all codes less than 50 are rejected 
Since they contribute no appreciable amount of precipitation. 
The codes used had values ranging from 50 to 99. Each of 
these codes is associated with a term made up of some combi- 
nation of three rainfall parameters. These are x, associated 
with light, continuous rainfall; y, associated with moderate, 
continuous rainfail; and z, which is a measure of heavy, 
continuous rain. These three parameters satisfy the follow- 


mid criteria: 


y = nx mm where l<n<5 
Zz = my mm where l<m<3 


z = 10 mm (in 3 hr) 


For this study x was set equal to 1.85, y = 5.66, and 
z = 8.13, the same values as derived by Tucker from his 
meeguency distributions of the codes recorded and the actual 
fod)! at Jjand stations. Thus, the amount of rainfall, 
in millimeters, associated with each current weather code 
can be determined from one of the above, or a combination or 
ratio of these terms. For example, weather code 64, which 
is "rain, not freezing, intermittent, but heavy at the time 
of observation", is described by the term z/2. Likewise, 
code 84, "showers of rain and snow mixed, moderate to heavy’, 


1S represented by the formula, ve : 


a2 






For the purposes of this study, the precipitation was 
calculated every three hours and then summed to obtain the 
total precipitation for that particular day. This was then 
subracted from the evaporation term for that day to yield 


the E-P time series for 1968 through 1970. 


B.. COMPUTATION OF THE ENERGY SPECTRUM 

The covariance function, Real Si) of a time series of 
random data is a measure of the general dependence of the 
values of the data at one time, t, with the values at 
another time, ttr. The covariance equation in integral form 
nS : 


R (ct) = 


: 
. a | cae iad (14) 
0 


] 
T +0 
This yields the average of the product of the value of a 
data point, x(t), with that of x(tt+tt), a point at some 
lagged time increment, +t. AS T approaches infinity, the 
resulting average product approaches the exact auto- 
covariance function. 

Since this data forms a finite data set, a lag window 
is applied to the covariance function. The Parzen window 


used in this analysis is defined as: 


P(t) = 1-6(r/T_)* + 6(t/T_.)> + = 01,2 --- T/2 (15) 
: (1-(x/T,))° eat 72 4 1, -=-1, (16) 
= 0 eV 1 (17) 
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where +t is the lag and fe ismeneetotdislag time. his 
window has the desirabte effect of minimizing any numerical 
instability due to side lobe resonance (Bendat and Piersol, 
971). 

The energy-density spectrum is obtained by Fourier 
transforming the covariance function. The form of the 


energy-density spectrum, corrected by the Parzen window is: 


fee] 


Sef) = P(t) RL (x) a7 1élfty, (18) 


with f being a specific frequency. The energy-density 
spectrum is the statistical estimate of the energy-density 
distribution for the frequency bands forming the total 


Spectrum. 


i COMPUTATION OF THE CROSS-SPECTRAL DENSITY 

When comparing two-time series of the same length and 
mumemincerval, the cross-covariance function similar to the 
covariance function, is used. For a given lag time interval, 


t, it has the form: 


Xy if 


T+ 


; 
ecjeestin — { oy (ey ae (19) 
0 


Because the data sets are records of finite length, the 

Parzen lag window is applied. Thus, similar to the covariance 
function and the spectral density, the cross-spectral density 
is the Fourier transform of the Parzen window lagged, cross- 


covariance function. 
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The cross-spectral density is a non-zero function, 


having real and imaginary parts of the form: 


L2) 


Syy (Ff) = | P(r) Sayed e léntt dt (20) 


= 00 


which can be expressed in terms of its real and imaginary 


parts: 


Syy(f) = Cyy(f) - id, (F) (21) 


where CY (f), the real co-spectrum, is given by: 


y 


yh) = 2 i LS ay ee - Scop SO COS I 2aey Molt 22) 
0 


and the imaginary quadrature spectrum, Q.. (f): 


Qy(f) = 2 { [Syy(t) = Syy(-r)] sin (2nfr) dr (23) 
0 


G. COMPUTATION OF THE COHERENCE FUNCTION 
The coherence function, for two records x and y, is a 
measure of the statistical dependence of two series. The 


coherence is defined as: 


1? ig [s. (Ff) 1° = 
xy ig it vy f 
where: 
ieee sn) <1 
XY 
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If vey (F) = 1, for all frequencies, then the two records 

are totally coherent and statistically dependent. If 

¥° (f) = 0 for all frequencies, then the x and y records 

wl totally independent. For intermediate range coherences 

three conclusions may be inferred concerning the two records 
(Bendat and Piersol, 1966); (1) the systems are non-linear, 

(2) noise is present in the measurements, or (3) the output 

response is due to more than one input function. It should 


be noted that this is the coherence squared, but will here- 


after be referred to as the coherence. 


ile COMPUTATION OF THE PHASE ANGLE 

The average angular difference by which the cross- 
correlated points in y(t) lead their corresponding components 
of x(t) is called the cross-spectral phase angle, Song the 
It is computed from the co- and quad-spectra by: 


(f) = arctan Boe / Sap eS), 


© xy 
Since this study deals with a time step of one day, the 
phase angles can be converted to a lag in days from: 

ey epee 0401745 
days pS SS (26) 
Parantt 

When either the value of the phase angle (measured in 
degrees) or the lag time (in days) is positive, then y(t) 
lags behind x(t) by that amount; if the value is negative, 
then x(t) lags behind y(t). 
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In the computations involved in this study there were 
8766 data points. A standard lag of 876, 10% of the total 
time series,was used. The 90% confidence interval for the 


energy densities was computed by the formula: 


when N = 8716, M = 876. 


The limits were found from tables to be: 


a 24.075 


b 5) SNe 


where a and b are the values which when subtracted and 
added, respectively, to the value of the spectral energy 


density result in the lower and upper 90% confidence limits. 
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A. REGRESSION ANALYSIS 

The regression analysis phase of this study was limited 
to the three-year time series of the NODC hydrocast data due 
to the desire to utilize the salinity data present in this 
data set alone. The period covered was from 1968 through 
1970, the same time span of data analyzed by Hansen (1973). 
At this point it is appropriate to visually inspect the 
Surface temperature and salinity data sets in a determinis- 
mie approach. 

Figure 4 shows the monthly means of the sea-surface 
temperature for 1968 through 1970 compared to Robinson's 
(1971) twenty-year long-term monthly means. A characteristic 
annual cycle is evident in all three years as well as the 
long-term mean. The only major difference from the long- 
term mean is a persistent cold anomaly which is present in 
the spring of all three years. 

The year 1968 was the most anomalous, being both colder 
in March and April and warmer for the balance of the year. 
Namias (1971) attributed this to bulk heating, advection 
and mixing. The year 1969 was more normal, except for an 
anomalous cooling in the summer, attributed by Roden (1970) 
to a southward displacement of the transition zone. This 
has been supported by Wickett and Thompson (1971a, 1971b) 


whose meridional transport data, Figure 5, shows a transport 
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of cooler water equatorward during April through October, 
1969. Except for the spring months, 1970 was more typical 
to the long-term mean than the preceding two years. 

Figure 6 shows the monthly means of surface salinity 
for 1968 through 1970. The most obvious fact from this 
figure is that the salinity shows evidence of having a 
semi-annual periodicity. In addition to this, there is the 
feature in 1970 of a reduction in the salinity, possibly due 
to the advection of abnormally fresh, cool subarctic water, 
indicating an equatorward shift in the subarctic transition 
zone. 

As stated previously, Hansen (1973) did much work on 
the same data series. Figure 7 is the result of his 
bivariate analysis of the monthly mean surface temperatures 
and salinities. It shows that the temperature follows four 
phases, namely the seasonal phases of summer, autumn, winter 
and spring, with autumn and spring appearing as transitional 
phases. Hansen also performed a regression test of the 
seasonal surface temperature versus the surface salinity and 
also the temperature at 50 m. The correlation coefficients 
are given in Table I. For most cases the correlations were 
poor, especially between the surface salinity and the surface 
temperature. The only good correlations occurred during the 
transition periods when both the surface temperature and 


Salinity were responding to the changes of the annual cycle. 
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Figure 7. Monthly Temperature and Salinity Modal Points 
Resulting from the Bivariate Analysis of the 
Data, 1968 through 1970 (from Hansen, 1973). 
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Table I 


Correlation coefficients, by season, resulting from a linear 
regression analysis of surface temperature, surface salinity, 
and 50 meter temperature data, 1968 through 1970 (from 
Hansen, 1973). 

1968 1969 1970 Mean 


She S00 m “Ste 50 m Sic 50 m SFE 50m 
Sa] Temp Sal Temp Sal Temp Sal Temp 


Winter UeOle O750) Ono ll 0.56 Oey, Oc Re miomtO 2.5 2 
Warming Uo wee 90 Oey 7) 0.70 8-02.11 Cay Gace 0 5/'5 
Summer S0eomeo. Sle 0e16 0.40  =0705 -0.29 -0.05 0.14 
Cooling Upovmeorvco 0204 0.98 -0701 OS SKS 0.47 0.94 


Since the surface salinity and temperature were 
obviously responding to different causal parameters, a 
further regression analysis was performed. The surface 
Salinity was regressed against many of the classical para- 
meters associated with its fluctuations. In this instance, 
there were seven variables chosen; the vapor pressure 
difference, E-P, the heat of evaporation, the wind speed, 
evaporation, precipitation and surface pressure. No separa- 
tion into seasons was made, all computations being made as a 
three-year series of daily average observations. 

The results of the regression study are shown in Table Ii 
Two major observations are apparent from this study. First, 
the standard parameters often associated with a change in 
Salinity, especially E-P, showed very little correlation with 


the surface salinity. However, this is quite reasonable 
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Table II 


Correlation coefficients, all parameters vs surface salinity. 


PARAMETER 1968 1969 1970 
ee. 094 043 gO F7. 
E-P 045 =o) -.016 
QE 02S on -.00] 
W -.030 .163 -.043 
ES O30 S03 TOUS 
B 034 . 166 004 
p -.030 -.019 .028 


considering the location of OWS NOVEMBER with its year round 
excess of evaporation over precipitation. The other major 
observation is that only one parameter is consistently 
highly correlated, i.e., the surface pressure. 

This latter fact could mean that the changes in the 
Surface salinity are largely due to an advective process, 
as 1S Supported by the fact that the other non-conservative 
terms associated with the air-sea interface have such poor 
correlation coefficients. This suggests that the changes in 
the surface salinity are related to some geostrophic process, 
most likely an advection of water from the subarctic region, 


corresponding to a southward shift of the transition zone. 
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B. SPECTRAL ANALYSIS 

The major thrust of this study was to use spectral 
analysis aS a means Of numerically describing the inter- 
actions at the air-sea interface at OWS NOVEMBER. Specifi- 
cally, it was to identify those meteorological parameters 
most closely associated with fluctuations in the sea-surface 
temperature. Since this latter term has a dominant annual 
cycle and most of the other parameters were compared with 
this term, the results presented here will primarily be 
those for the yearly cycle. The logic of this is further 
demonstrated by comparing the percentages of the variance 
explained by the yearly cycle for the different terms. The 
percentages shown in Table III were calculated from the 
energy-density functions which were produced for each time 
series. This function is a measure of the variance contained 
at that frequency. To compute the variance contributed Dy 
the yearly cycle a bandwidth was specified by those frequencies 
for which the value of the energy-density function decreased 
to less than half the value at the yearly cycle. 

From Table III it can be seen that 86% of the total 
variance of the sea-surface temperature is explained by the 
yearly cycle alone. For this reason, even though the 
Spectral results were computed for frequencies as short as 
the monthly cycle, the emphasis of this study was placed on 
the yearly cycle. Also, since the major areas of interest 


were the correlations with the heat budget terms, each of 
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these.is discussed describing its influence on the sea- 
Surface temperature. In addition, the contributing effects 
of the parameters involved in the computation of that heat 


flux term will be discussed. 


Table III 


Contribution of the yearly cycle to the total variance explained 


Parameter % Of Total Variance 
SS] 86 
Ie 80 
lige - Us 20 
Qe, 86 
e. 46 
ene ca 14 
E (evaporation) 34 
W 8 
we 9 
Cl (cloud Lovie 7 
S 8 
QS 70 
QB 87 
Gir 11 
QC V7 
QN 42 
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The first heat budget term to be examined is the solar 
radiation, QS, which like the surface temperature exhibits 
a strong annual periodicity. From Tabe III, it can be seen 
that 70% of the total variance of this term is explained by 
the annual cycle. The term which primarily controls the 
variation of QS is the cloud cover, Cl, and in Table IV the 
Spectral results for the yearly cycle are given for the 
miceractions of the SST and Cl with QS. Along with a 
correlation of the cloud cover with the SST, the complete 
spectral results are shown in Figures 8, 9, and 10. Figures 
A-1, A-2, and A-3 of Appendix A show the respective auto- 
correlation functions and Figures B-1] and B-2 of Appendix B 
illustrate 15-year segments of the 24-year time series of 
the SST and QS. 

As would be expected, Table IV shows that the solar 
radiation is strongly coherent with the SST at the yearly 
cycle. However, there is a lag of almost 60 days, which 
cannot be explained by the cloud cover correlation with its 
low coherence and therefore dubious lag time. Also, there 
is a strong coherence at the semi-annual cycle, which appears 


to reflect the effect of the cloud cover as seen in Figure 10. 


ade slay 


Spectral results for the solar radiation, 
SST and cloud cover at the yearly cycle. 


Parameter Coherence Phase Lag 
Codes (degrees) (days) 

mS vs SST ~94567 - 60.0 -58.5 

WS” V Sane | Boo 230 aO.. 2 We2 510 
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Figure 8. Spectral density, coherence and phase for 
Eiemsonevs the solar radiation. 
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Figure 9. Spectral density, coherence and phase for 
the SST vs cloud cover. 
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This rather large lag is equally evident in Table V 
which shows the months tn which SST and the heat transfer 
terms achieve their maximum and minimum values. Along with 
Figures B-] and B-2, it can be seen that the surface 
temperature is a maximum during September and at a minimum 
during March. However, the solar radiation peaks in August 
and iS a minimum in January. This is approximately the 58 
day lag shown in Table IV. This lag can be explained by the 
fact that the solar radiation is absorbed and transferred 
throughout the upper mixed layers of the water column. Thus, 
a lag of two months exists before the total effects of the 
solar radiation are represented by a warming or cooling of 


the surface temperature. 


Table V 


Months of maximum and minimum value for the 
surface temperature and the heat budget term. 


Parameter Minimum Maximum 
SST March September 
QS | January August/September 
OE May November 
QB July/January March/November 
QC August March 
QN December June 
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In Figures 8, 9, and 10 there are also noticeable peaks 
at the six month cycle and three month cycle. These results 
are thought to be, for the most part, artifacts of the yearly 
cycle and are relatively small when compared to the yearly 
value; yet, as mentioned before, several show evidence of 
being real. 

The major result of this portion of the study is that 
miesuse of the new formula for QS results in Significantly, 
higher estimates of this term than previously calculated. 
Fluctuations in solar radiation precede those of the surface 
temperature by approximately two months and thus the solar 
radiation is a good predictor of the surface temperature, 
as would be expected. 

The back radiation, QB, represents the loss of heat 
energy as long-wave radiation. It is strongly dependent 
iagaeeene SuYtace temperature and the vapor pressure of the 


atmosphere. Table VI summarizes the results for the yearly 


cycle of QB and associated terms with the surface temperature. 


Table VI 


Spectral results for the back radiation and 
related parameters at the yearly cycle. 


HVS SST 5 EAS leo ES 18 
is ve) SST .98800 - 6.0 sie gre, 
a ViSae So) .67194 -49.] -47.9 
7 VS S51 + eee is) - 6.4 - 6.2 
QB vs Cl 77005 Dome. = On| 
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From Table III it was found that about 87% of the total 
variance of the back radiation was contained in the yearly 
cycle; yet, Table VI shows only a slight coherence between 
the SST and QB for the yearly cycle and a large dubious lag 
value. This can be explained by examining the interactions 
involved. The back radiation term decreases as the SST 
increases due to the effect of the humidity in the atmosphere 
directly above the surface. Thus, over the yearly cycle, 

QB would ideally be 180 degrees out of phase with the surface 
temperature. However, the back radiation is also decreased 
by an increase in the cloud coverage. Thus, the condition 
illustrated by Figure B-3 and Table V may be seen where QB 

is shown to have two maxima and minima on roughly a semi- 
annual basis. From Figure A-4 it may be seen that this 
feature is superimposed over a larger annual cycle which 
explains the high percentage of the variance contained by the 
yearly cycle. 

The complete spectral results for the above correlations 
are given in Figures 11] through 15. From Figure 11 it is 
seen that QB has two peaks, one at the yearly cycle and 
another lesser one at one-half the yearly cycle. Of the 
related parameters, only the air temperature and the vapor 
pressure of the air are highly correlated and then at only 
the yearly cycle. The air-sea temperature difference is only 
mildly correlated at the yearly cycle and its associated lag 


is highly suspect for that reason. The air-sea temperature 
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Figure 11. Spectral density, coherence and phase for 
mew ootevs the back radiation. 
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Figure 12. Spectral density, coherence and phase for 
the SST vs air temperature. 
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Figure 13. Spectral density, coherence and phase for 
the SST vs the air-sea temperature difference. 
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Figure 14. Spectral density, coherence and phase for the 


SSievocetie vapor pressttre of the air. 
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Figure 15. Spectral density, coherence and phase for the 
back radiation vs cloud cover. 
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difference demonstrates a coherence at the 70-day cycle. 
Also the back radiation and cloud cover show evidence of 
this possible seasonal coherence. 

The major conclusion to be drawn from this part of the 
study is that the back radiation is semi-annual with an 
underlying yearly cycle. As a predictor for the surface 
temperature, the back radiation is poorly associated with 
meuctuations tn SST. 

The next heat budget term examined was the heat of 
evaporation, QE. Only 11% of the total variance of QE is 
explained by the yearly cycle which in itself means that QE 
is a poor predictor of the surface temperature. This is 
verified by Table VII which shows only a moderate coherence 
for QE and the SST at the yearly cycie. This iessens the 
credibility of the large lag shown; also from Table V it can 
be seen that the actual lag is more of the nature of two 
months. Additionally the lag is positive which means that 


fluctuations of QE follow those of the SST. 


Table VII 


Spectral results for the heat of evaporation 
and related parameters at the yearly cycle. 


Parameter Coherence Phase Lag 
Codes (degrees) (days) 
ME VS SST .84783 SiiD) 95.0 
ee, vs Sst -86425 aes A Tara 
evs SST . 9559 3 roe. 135.4 
Wo vs SST 95512 135.4 132.0 
Eevs. SST ae 76.9 an 
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Whereas QE is only weakly correlated with the SST at 
the annual cycle, the wind speed and wind speed squared are 
highly correlated to the SST with large positive lags 
(approximately 4 months). These would seem to explain much 
of the lag in QE when they are coupled with the vapor 


pressure difference. QE is a transfer mechanism of the form: 


=b 


maven. 2 


A (28) 


Q 


Z 


where: 
a is the eddy coefficient for water vapor 
at is the humidity gradient above the sea surface. 


Both of these are highly dependent upon the wind speed. 
Even though the vapor pressure aifféerence and evaporation 
are weakly correlated to the SST, they combine with the wind 
speed to produce the low coherence and large lag for QE. 
Figures 16 through 20 give the complete spectral results 
for the correlations shown. Being computed from the surface 
temperature, e and SST are naturally perfectly correlated 
and is not included. The yearly cycle dominates in all 
except for QE, the vapor pressure difference and the evapora- 
tion, as would be expected. Both the wind speed and wind 
Speed squared show a definite coherence at about the 80-day 
cycle. The autocorrelation functions are found in Figures 


A-9 through A-13 and a 15-year plot of QE and its anomaly 


1S shown in Figure B-4. 
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Figure 16. Spectral density, coherence and phase for the 


SST vs the heat of evaporation. 
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Figure 18. Spectral density, coherence and phase for the 
SST vs the wind speed. 
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Figure 19. Spectral density, coherence and phase for the 
SST vs the wind speed squared. 
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Figure 20. Spectral density, coherence and phase for the 


SST vs the evaporation. 
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The heat of conduction, QC, is the smallest of the heat 
flux terms and only 17% of its total variance is explained 
by the yearly cycle. From Table VIII it can be seen that 
this term also has an unusually large lag, approximately 
130 days, and the lowest coherence of any of the heat budget 
terms with the surface temperature. From Figure 21 it can 
be seen that QC is nearly of the form of white noise with 
no distinct periodicities. This is confirmed by Figure B-5 
which shows that QC is constantly small and contains few 
anomalies of note. 


Like QE, QC is a flux term of the form 


oT 
QC = sla ae (29) 
wnere 
An is the eddy coefficient of heat conduction, and 
at is the temperature gradient of the air overlying 


the sea surface. 


For an area of the nature of OWS NOVEMBER, dominated by 
fairly constant temperature patterns for both the air and 
the ocean, it can be seen that this term would be expected 
to be fairly constant. It does contain quite a deal of 
energy in cycles other than the annual and shows evidence 
of good coherence at the semi-annual, seasonal and monthly 


periodicities. 
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Figure 21. Spectral density, coherence and phase for the 


SST vs the convective heat. 
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TABLE VIII. 


Spectral results for the convective 
heat vs the SST at the yearly cycle 


Parameter Coherence Phase Lag 
Codes (degrees) (days) 
OC iiss SST . 75694 134.1 130.7 


The final heat budget term, the net heat transfer, QN, 
is perhaps the most interesting of the heat equation terms, 
representing the gain or loss of heat by the water column. 

To this point, the advective and local heat exchange terms 
have been ignored. However, since the values for the net 
radiation computed from Equation (1) did not result in zero 
values, then it is evident that these terms have to be 
considered. If QN is then viewed as a combination of the 
advective heat term and the local generation of heat it takes 


on the form: 


QN = Q) + Q, 


where: 
oe 1s the advection of heat by currents, and 


Q» is the local heat generated. 


Over the annual cycle Q, = 0; therefore, QN must be 
representative of only the advection term. This is Supported 
by the fact that for the entire 24-year period QN was found 
to have a slope of only 0.003 which signifies that no long 


term generation of heat has occurred. 
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From Tables III and IX it can be seen that the net heat 
exchange is dominated by the annual cycle and is highly 
coherent with the surface temperature at this period. 
Forty-two percent of the total variance is explained by the 
yearly cycle, which allows for some seasonal effects and | 
Figure 22 shows some evidence of high coherence at the 
semi-annual and seasonal cycles. The lag is confirmed by 
both Table V and Figure B-6. Much of this lag is undoubtedly 
due to the effect of the 58-day lag of the solar radiation 
and the balance by the other heat budget terms. In addition, 
the lag of QN and the surface pressure, considered here 
because of the previously noted high correlation of pressure 
with the surface salinity, may also contribute to the lag of 


QN and SST. 


TABLE IX 


Spectral results for the net radiation 
and related parameters at the yearly cycle 


Parameter Coherence Phase Lag 
Codes Squared (degrees) (days) 
QN vs SST ~97954 ae Lon - 74.8 
ZoevS SO! . 79993 -120.6 -118.6 
ON vs PS ee 0 - 46.9 - 45,3 
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Spectral density, coherence and phase for the 
SST vs the net radiation. 
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Figures 22 through 24 give the-.complete spectral 
results for the net heat exchange, the surface temperature, 
and the surface pressure correlated versus each other. In 
all three instances the yearly peak dominates. 

The major result of this segment is that the new solar 
radiation formula results in larger positive values of QN. 
Secondly, for the yearly cycle QN is a measure of heat 
advection and at OWS NOVEMBER this advection is of nearly 
equal importance as the solar radiation as an influence on 
the surface temperature. The negative lag makes QN a fairly 
good predictor of the surface temperature, even though much 
of the variance of the net heat exchange is contained in the 
seasonal cycles. The surface pressure seems to be less 
correlated with the surface temperature than it is with the 
Salrinity. 

Very little has been mentioned about the autocorrela- 
tions; however, White (1974) has used this function in the 
Study of SST and baroclinic transport anomalies.. He found 
an apparent 5-year periodicity in the SST anomalies in the 
Meinity of OWS NOVEMBER. This cycle was attributed to 


deeper anomalous features in the baroclinic structure. 
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Figure 24. Spectral density, coherence and phase for the 
Netmicadioarilom VS the surface pressure. 
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iy. CONGIUS LONS 


The aim of this study was to use the tools of spectral 
analysis and other statistical approaches to describe the 
interactions of the air-sea interface at OWS NOVEMBER, 
specifically as to finding those parameters most closely 
associated with fluctuations of the daily average surface 
temperature and surface salinity. 

The results of this research showed that: 

il. The surface salinity from 1968 through 1970 was 
shown through regression analysis to be poorly correlated 
with the parameters traditionally associated with changes 
in the observed salinity. The evaporation, evaporation 
minus precipitation, and precipitation were found to have 
little correlation with fluctuations in surface salinity. 
These terms are non-conservative terms associated with the 
air-sea interface and their low correlation coefficients 
would appear to Racce chatr importance in controlling the 
mr ace Salinity at this location in the Pacific. 

Zz. The surface pressure was found to be the most 
strongly correlated of al! the parameters regressed against 
the surface salinity. It was therefore hypothesized that 
this represents the presence of a movement of a baroclinic 
system of cooler water from the subarctic region into the 
Vicinity of OWS NOVEMBER, corresponding to a shift southward 


memcne transition zone. 


iS 





S. For the yearly cycle, it was shown that the net 
radiation is a measure of the advection present and that 
this advection is of nearly equal magnitude as the solar 
maa1ation. 

4, The new formula for the solar radiation leads to 
values for this term approximately 20 to 30 percent higher 
than previously computed. In turn, this leads to much 
larger positive values of the net heat exchange term and 
larger positive anomalies for both terms. 

or From the spectral analyses it was found that all 
the 24-year time series were dominated by the yearly cycle; 
except for the back radiation, cloud cover, wind speed and 
wind Piwed Squared series which showed evidence of semi- 
annual and seasonal cycles as well. 

6. As predictors of the surface temperature, the 
Solar radiation and advection in the form of the net heat 
exchange term were the only terms of suitable phase and 
coherence at the yearly cycle to be valid. The other heat 
transfer terms seems to have less effect upon fluctuations 
of the surface temperature. Several of the observed para- 
meters, such as the air temperature, were also highly 
coherent and, on the whole, had smaller lags associated with 


them. 
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APPENDIX A 
THE AUTOCORRELATION FUNCTIONS 
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Figure A-]. 
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Autocorrelation curves for the SST and the solar 
radiation, QS (abscissa is given in hundreds 
of days). | 
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Figure A-2. Autocorrelation curves for the SST and cloud 
coverage (abscissa is given in hundreds of days). 
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Figure A-4. Autocorrelation curves for the SST and the back 
radiation (abscissa is given in hundreds of days). 
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Figure A-5. 
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Autocorrelation curves for the SST and the air 
temperature (abscissa is given in hundreds of days). 
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Figure A-6. Autocorrelation curves for the Soll and the 
air-sea temperature difference falbiste ic sia 119 


given in hundreds of days). 


83 





Q2s 


Dee 


Wa) 

a 

2 | 

fom] 

of 

[at J 

4 

:? 

fe), 

Gr 

v? |! 

(ea ‘ aan 

oc \ ia: \ nm ay 034 
KK e@ 
a 

44 

cer: 

a 

| wes ST 
c? 

=| 

o | 

4 


madre A-7. Autocorrelation curves for the SST and the vapor 
pressure of the air (abscissa is given in 
hundreds of days). 
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Figure A-8. Autocorrelation curves for the back radiation 
and cloud cover (abscissa is given in hundreds 


of days). 
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Figure A-9. Autocorrelation curves for the SST and the heat 
of evaporation (abscissa is given in hundreds 


of days). 
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Figure A-10. Autocorrelation curves for the SST and the vapor 
pressure difference (abscissa is given in 
hundreds of days). 
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Figure A-11. Autocorrelation curves for the SST and the 
wind speed (abscissa is given in hundreds 
Siedays ). 
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Figure A-12. Autocorrelation curves for the SST and the 
wind speed squared (abscissa is given in 
hundreds of days). 
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Figure A-13. 
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Autocorrelation curves for the SST and 
evaporation (abscissa is given in hundreds 
of days). 
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Figure A-14. Autocorrelation curves for the SST and the 
convective heat (abscissa is given in hundreds 
Oneedays ). 
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mrcuire A-1t5. Autocorrelation curves for the SST and the net 
radiation (abscissa is given in hundreds of 
days). 
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Figure A-16. Autocorrelation curves for the SST and surface 
pressure (abscissa is given in hundreds of days). 
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Figure A-17. Autocorrelation curves for the net radiation 
and the surface pressure (abscissa is given 
in hundreds of days). 
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- APPENDIX B 
PIFIEENSrEARK TIME SERTES SEGMENTS 


Fifteen-year time series of the 
surface temperature and its anomaly 


Fifteen-year time series of the 
solar radiation and its anomaly ---- 


Fifteen-year time series of the 
back radiation and its anomaly ----- 


Fifteen-year time series of the 
heat of evaporation and its anomaly 


Fifteen-year time series of the 
convective heat and its anomaly ---- 


Fifteen-year time series of the 
net radiation and its anomaly ------ 
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